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Abstract 



Wc report new features of the typical mixed-morphology (MM) supernova remnant (SNR) W44. In the 
X-ray spectra obtained with Suzaku, radiative recombination continua (RRCs) of highly ionized atoms are 
detected for the first time. The spectra are well reproduced by a thermal plasma in a recombining phase. 
O-f The best-fit parameters suggest that the electron temperature of the shock-heated matters cooled down 

q . rapidly from ~ IkeV to ~ 0.5 keV, possibly due to adiabatic expansion (rarefaction) occurred ~ 20,000 

years ago. We also discover hard X-ray emission which shows an arc-like structure spatially-correlated 
with a radio continuum filament. The surface brightness distribution shows a clear anti-correlation with 
12 CO (J = 2-1) emission from a molecular cloud observed with NANTEN2. While the hard X-ray is most 
likely due to a synchrotron enhancement in the vicinity of the cloud, no current model can quantitatively 
predict the observed flux. 

' Key words: ISM: abundances — ISM: individual (W44) — supernova remnants — X-rays: ISM 

CO" 
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1. Introduction (Ohnishi ct al. 2011) and W28 (Sawada & Koyama 

2012). The strong RRC emission can be observable when 

qq ■ Mixed- morphology (MM) supernova remnants (SNRs) an ionization temperature (feT^ 1 ) is significantly higher 

are classified based on the characteristic morphology of than electron temperature (kT c ) so that the free-bound 

CN1 ■ centrally-filled thermal X-ray emission with a synchrotron transition (recombination process) becomes dominant — 

radio shell (Rho & Petrc 1998; Lazendic et al. 2006; Vink hereafter, we call it a recombining plasma (RP). These 

2012). More than 25% of the X-ray-detected Galactic discoveries are dramatically changing our understanding 

SNRs are categorized as the MM type (Jones et al. 1998). of the SNR evolution, since most young or middle-aged 

Most of the MM SNRs indicate evidence for interaction SNRs have either an ionizing plasma (IP: kT c > kT z ) or a 

with dense interstellar medium (ISM), CO and Hi emis- collisional ionization equilibrium plasma (CIE: kT z = kT e ). 

sions or OH (1720 MHz) masers in some cases (Rho & While hints of RPs were reported from a few non-MM 

Petre 1998). They are also associated with TeV/GeV 7- SNRs (Hughes et al. 2003; Broersen et al. 2011), all the 

ray sources, suggesting the presence of molecular clouds robust results have been found only in MM SNRs. It 

in their vicinity (e.g., Albert et al. 2007; Aharonian ct al. is suggested, therefore, that the formation process of the 

2008; Hewitt et al. 2009; Abdo et al. 2010). Such ambient RPs is closely related to the characteristic environment of 

environment would play some roles in the evolutions of the the MM SNRs. 

morphology and plasma condition in these SNRs; e.g., the The environment may also affect the particle accelera- 

central X-ray emission might be enhanced by cool clouds tion in SNR blast waves. The shock waves propagating 

left relatively intact after the passage of blast waves to into the inhomogeneous ambient medium may produce 

slowly evaporate in the hot SNR interior (White & Long strong turbulence in the magnetic field, which enhances 

1991), or density gradients in the pre-existing ISM could the efficiency of the particle acceleration in this region 

contribute the emission profile (Petruk et al. 2001). (e-g-: Giacalone & Jokipii 2007; Inoue et al. 2012). 

Recent Suzaku observations revealed the presence of en- In this paper, we focus on one of the most typical MM 

hanced radiative recombination continua (RRCs) in X- : 

r 1 n «-n 1 n»Tn t/-i a An /-%r • ■ kT z represents populations of ionization states. In collisional 

ray spectra of several MM SNRs, IC443 (Yamaguchi ionization equilibrium, fcT z is equal to the plasma temperature; 

ct al. 2009), W49B (Ozawa et al. 2009), G 359.1-0.5 see also Masai (1994). 
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Table 1. 


Observation logs. 








Name 


Obs. ID 


Obs. Date 


(R.A., Dec.) J2000 


(l,b) 


Exposure 


Source 
Background 


W44 
GR 


505004010 
500009020 


2010-Apr-10 
2007-Aug-18 


(18 h 56 m 08 s , +01°23'19") 
(18 h 44 m 01 s , +04°04'41") 


(34.°7, -0.°41) 
(28.°5, +0.°21) 


61.1 ks 
98.9 ks 



SNRs, W44 (also known as G 34.7-0.4 or 3C 392) to study 
its X-ray properties and their relation with the environ- 
ment. W44 is a middle-aged (~ 20, 000 yr; Wolszczan et al. 
1991) SNR located on the Galactic plane at a distance of 
D ~ 3kpc (Claussen et al. 1997). In the radio band, W44 
shows a distorted shell of the synchrotron emission (e.g., 
Kundu & Velusamy 1972; Handa et al. 1987). Recent 
high-resolution radio observations revealed more compli- 
cated filaments in the radio continuum shell (Giacani et al. 
1997; Castelletti et al. 2007). Molecular radio emission 
was detected from the vicinity of W44 by several obser- 
vations of the CO (Wootten 1977; DeNoyer & Frerking 
1979; Seta et al. 1998; Reach et al. 2005) and Hi (Koo 
& Heiles 1995) lines. More recently, Seta et al. (2004) 
showed a whole view of the high-resolution 12 CO emis- 
sion map with the 17"-beam of the Nobeyama 45-m radio 
telescope and reported many hints for the shock-cloud in- 
teraction, including cloud evaporation. The presence of 
1720 MHz OH maser emission is also a sign of the interac- 
tion between the SNR shocks and the ambient molecular 
clouds (e.g., Claussen et al. 1997; Frail & Mitchell 1998). 

In the X-ray band, W44 was found to have center- 
filled morphology by ROSAT observations (Rho et al. 
1994). Based on the Chandra result, Shelton et al. (2004) 
proposed that the centrally bright X-ray morphology is 
caused by either entropy mixing in the form of ther- 
mal conductions, or evaporation of swept-up clouds. The 
analyses of combined X-ray spectra of ROSAT, EXOSAT 
(ME), Einstein (SSS), and Ginga revealed that the plasma 
was IP and has not reached CIE, and the spectra are in- 
terpreted by cloud evaporation (Rho et al. 1994; Harrus 
et al. 1997). ASCA obtained higher-quality spectra and 
determined the ionization temperature (kT z ) of Si from 
the flux ratio of Ka lines between H-like and He-like ions 
(H-like/He-like). The result was that kT z is nearly equal 
to the electron temperature kT e (Kawasaki et al. 2005), 
or the thermal plasma was nearly in CIE. 

This paper reports revised X-ray features of W44 made 
with the Suzaku observation. The bright X-ray flux and 
apparent X-ray size of ~ 0.°6 in diameter are best suited 
to study the detailed plasma structure. The high sensitiv- 
ity and the good energy resolution of Suzaku enable us to 
determine the spatial structure of the X-rays accurately, 
in particular, that of the thermal plasma. We report the 
first detection of an RP in this SNR and discovery of hard 
X-ray emissions. The evolution of the RP is quantitatively 
examined, and then some scenarios for its origin are ad- 
dressed. The hard X-ray structure is compared to that 
of the molecular cloud with the NANTEN2 observations. 
The origin is discussed in the context of a shock-cloud 
interaction. Throughout this paper, the distance to the 
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Fig. 1. Top Left: Three-color XIS image of W44 overlaid 
on the Galactic coordinate. The red, green and, blue corre- 
spond to the 1.0-2.0 keV, 2.0-5.0 keV and 5.0-7.0 keV bands, 
respectively. The white contours represent X-ray surface 
brightness obtained by the ASCA GIS observation (Kawasaki 
et al. 2005). The FoV of the XIS is shown with the white 
square. The white cross indicates the location of the pulsar 
PSRB1853+01. Top Right: Hard X-ray band (5.0-7.0 keV) 
image obtained with the Suzaku XIS. The white contours rep- 
resent the 1442.5 MHz radio continuum emission with VLA. 
Bottom: Example of the source spectrum extracted from the 
region C. The solid line is the CXB model spectrum. 

SNR is assumed to be 3kpc, and the errors quoted are at 
the 90% confidence level. 

2. Observations and Data Reduction 

We performed a one-pointing observation of W44 with 
X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) 
on board Suzaku. The field of view (FoV) covers the X- 
ray brightest region in this SNR. In the southwest part 
of the FoV, the shock-cloud interaction has so far been 
confirmed by the CO observations (e.g., Seta et al. 2004). 
The observation log is presented in table 1. 

The Suzaku observation was performed in April 2010, 



Suzaku X-ray Observation of W44 



3 



during the A05 observing cycle. We used data of XIS1 
(back-illuminated CCD; BI CCD), and XISO and XIS3 
(front-illuminated CCDs; FI CCDs). XIS2 was unavail- 
able since 2006 November 9, possibly due to a microme- 
teorite damage. We used the revision 2.4 of the cleaned 
event data and combined the 3x3 and 5x5 pixel events. 
The calibration database (CALDB) updated in September 
2011 was used for the data reprocessing. We performed 
data reduction with the version 6.11 of the HEAsoft tools 
(version 18 of the Suzaku Software). 

We also performed observations of the 12 CO ( J = 2 
1) emission from 2010 August to 2011 January with 
NANTEN2, the 4 m-submillimeter telescope of Nagoya 
University at Atacama (4865 m above the sea level) in 
Chile. The half-power beam width at the frequency of 
the 12 CO (J = 2-1), 230.580GHz, was - 90". The ob- 
servations were made by on-the-fly (OTF) mapping mode 
and each spatial region was mapped several times in dif- 
ferent scanning directions to reduce scanning effects. The 
final pixel size of the gridded data was 30". The spectrom- 
eter was a digital Fourier transform spectrometer (DFS) 
with 16384 channels, providing a velocity resolution of 
0.08 km s -1 . The root mean square noise per channel af- 
ter smoothing was ~ 0.3 K. 

3. Analysis and Results 

In the following spectral analysis, we used the SPEX 
software version 2.02.02 (Kaastra et al. 1996). We em- 
ployed xisrmfgen and xissimarfgen to generate the re- 
distribution matrix files (RMFs) and ancillary response 
files (ARFs), respectively (Ishisaki et al. 2007). The en- 
ergy band around the neutral Si K-shcll edge (1.7-1.8 keV) 
is ignored, because an accurate response function in this 
band was unavailable at the time of the data analysis 2 . 

3.1. Background Subtraction 

The X-ray images after subtraction of the non X-ray 
background (NXB) (xisnxbgen; Tawa et al. 2008) are given 
in figure 1. The left and right panels show the three- 
color (three-energy bands) and hard X-ray images of W44, 
respectively. To make the source spectra, we further need 
to subtract the cosmic X-ray background (CXB) and the 
Galactic ridge X-ray emission (GRXE). 

The major component of the X-ray background is 
GRXE. Since the X-ray emission of W44 extends more 
widely than the FoV (figure 1), background data were 
not taken from the same FoV. We used the nearby sky 
on the Galactic ridge (hereafter GR). The observation log 
for the GR is also shown in table 1. The GR observa- 
tion was aimed 6.° 2 and 0.°2 apart from W44 along the 
Galactic longitude and latitude, respectively. This sepa- 
ration cannot be ignored for the background estimation, 
since the GRXE flux depends largely on the position of 
the Galactic plane, espetially on the Galactic latitude. 
The spatial distribution of the GRXE flux is modeled by 
a two-exponential function with scale lengths of 2500 pc 
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Fig. 2. The top panel shows the region A spectrum fitted 
with the NEIJ plus two-Gaussian model (see text). Black and 
gray represent the FI and BI data, respectively. The count 
rate of FI is multiplied by a factor of 10. The lower panels 
represent the residuals from the models of CIE (a), RP (b), 
NEIJ (c), and NEIJ plus two Gaussians (d). 

(17°) and 45 pc (0.°3) along the Galactic longitude and lat- 
itude, respectively (Uchiyama et al. 2011). Therefore, the 
GRXE flux around W44 is estimated to be ~ 40% of that 
at the GR position. We subtracted normalized GRXE 
spectrum, resulting no significant excess above 6.0 keV in 
the source spectra of the regions A and B. This is con- 
sistent with the imaging analysis (figure 1) where the no 
excess was found in the hard X-ray band, with an excep- 
tion of the region C. We hence conclude that the GRXE is 
properly subtracted. The hard X-ray excess in the region 
C is separately discussed later. 

The surface brightness of the CXB is less than ~ 5% of 
that of the GRXE. Therefore, possible spatial fluctuation 
in its flux can be negligible compared to the GRXE. We 
applied a power-low CXB model with a photon index of 
1.4, following the result by (Kushino et al. 2002). 

3.2. Fit with Thermal Plasma Models 

In figure 1-left, we see a bright centrally- filled structure 
in the soft X-ray band (1.0-5.0 keV). For spectral studies, 
we selected three representative regions, A, B, and C, in 
figure 1 left; the former two regions exhibit no hard X-ray 
emission, while the remaining one shows clear excess in 
the hard X-ray band. From the region A, the overlapped 
ellipse of the region C was excluded. 

We first analyze the spectrum of the region A. The 
previous observations suggested that the X-ray spec- 
trum of the SNR's center is characterized by a ther- 
mal plasma in either ionizing state or CIE (e.g., Harrus 
et al. 1997; Kawasaki et al. 2005). Therefore, we ap- 
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Fig. 3. Spectra of the regions A, B and C fitted with the NEIJ (+ two Gaussians) plus power-law model (see text). Black and gray 
correspond to the FI (xlO) and BI data, respectively. The residual is shown in each lower panel. 



Table 2. Best-fit parameters (see figure 3) 



Components 


Parameters 




Value 




Thermal Component (NEIJ) 
Abundance (solar) t 


kT cl [keV] 

kT c2 [keV] 

log(n e t [scm -3 ]) 

EM* [10 58 cm" 3 ] 

Ne 

Mg 

Si 

S 

Ar (=Ca) 
Fe (=Ni) 


region A 
1.00±0.05 
0.47±0.01 

11 7O+0.06 
iJ - 1 3 -0.03 

42.6 + ^ 3 

1 02 +oiM 

i - u *-0.Q5 

l-20±g;g| 

1.6818:85 

2 02+ - 10 

z - uz -0.16 
1 01 +0-34 

0.05±0.01 


region B 

0. 92±0.03 

n 49+0 05 

u -* a — 0.06 
n 77+0 .15 
LL - ' ' -0.29 

28.7+jj-J 

1 - za -0.32 

1 26+ 27 
n qe+o.28 

1. Ul_ 27 

<1.23 

0fi +005 


region C 

o.95l8;8I 

0.50±0.02 

11 oc+0.09 
ll.QO_ 25 

26. 2l 22 

i 97 +0.16 
L - z 1 -0.40 

1 46 +0 - 21 

1 - ao -0.13 

9 -| r+0.40 
z ' 1J -0.22 

2 59+ - 37 

z - dy -0.60 

<1.09 
<0.07 


Powerlaw Component 


r 

flux [lCT^crgscm^s- 1 ] 


2.12 (fixed) 
0.15 ±0.09 


2.12 (fixed) 
< 0.14 


2.12±0.11 

-i Q7+0.22 
-L" 3 ' -0.42 


Absorption 


7V H [10 22 cm" 2 ] 


1.3±0.2 


1.2±0.2 


1.3±0.2 




xVdof 


525/347 


323/219 


202/237 



* Volume emission measure at the distance of 3 kpc. 

' The other elements not listed here are fixed to solar values. 

The errors are in the range A% 2 < 2.7 on one parameter. 



plied a single-component CIE model with free parame- 
ters of plasma temperature kT c (= kT z ), emission measure 
(EM — J n c nnV), and column density Nn for interstellar 
absorption. Here, n , jih, and V are the number densi- 
ties of electrons and protons, and the emitting volume, 
respectively. The abundances relative to the solar values 
(Anders & Grevesse 1989) of Ne, Mg, Si, S, Ar and Fe 
were allowed to vary freely, while the Ca and Ni abun- 
dances were linked to those of Ar and Fe, respectively. 
Abundances of the other elements were fixed to one solar. 

The best-fit temperature of the CIE plasma model (kT c ) 
was obtained to be O.6OI8 02 keV, but this model was 
statistically unacceptable (x 2 /dof = 2121/353) with large 
residuals as shown in figure 2-a. The largest residual was 
found at 2.0 keV, which corresponds to the Si Lya emis- 
sion. Therefore, it is suggested that the ionization state of 
this element is higher than that in the best-fit CIE state. 



In addition, the hump- like were found around 2.7 keV and 
3.5 keV. These correspond to the K-shell binding energies 
of H-like Si (2.67 keV) and H-like S (3.48 keV). As is the 
case of IC443 (Yamaguchi ct al. 2009), these residuals are 
most likely to be RRCs in a recombining plasma (RP). 

We thus introduced a single-component RP model, 
where kT e and kT z /kT e were allowed to be free param- 
eters. Then the hump-like structures in the residual 
and the excess at the Si Lya were gone (see figure 2- 
b). The best-fit kT c and kT z /kT c were 0.42 ± 0.02 keV 
and 1.52 ±0.05, respectively However, the single RP 
model still failed to fit the overall spectrum with % 2 /dof 
of 1175/352, with a large discrepancy between the data 
and model in the low-energy band (< 2.0 keV). Among 
the residual features, a "dip" at ~ 1.5 keV (shown with 
an arrow in figure 2-b) is most remarkable. This corre- 
sponds to the Mg-Lya line, suggesting that the ionization 
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temperature kT z of Mg would be lower than those of Si 
and S. We supposed that the low-energy spectrum is dom- 
inated by a lower-temperature plasma attributable to the 
shocked ISM, and hence added a low-temperature CIE 
model. However this combined model (RP+CIE) gave no 
essential improvement in the low energy residuals. 

Finally, we consider that each heavy element has an in- 
dependent ionization temperature, as was first pointed out 
in IC443 (Yamaguchi et al. 2012). Sawada & Koyama 
(2012) also found that kT z for Ne and Mg in W28, an- 
other MM SNR, are closer to the CIE (kT z /kT e =1.2-1.7) 
than those for Si and S (kT z /kT e =2A-2.5). They success- 
fully demonstrated that this variation can be explained 
by the element-dependent recombination timescale (see 
also Smith & Hughes 2010). We therefore tried a non- 
equilibrium ionization jump (NEIJ) model for W44. This 
model describes a plasma state when the initial plasma 
is in CIE with the temperature of kT c i, then only the 
electron temperature dropped to kT e 2 by a rapid electron 
cooling in the past. Assuming a constant kT c2l the fol- 
lowing evolution of the RP is traced as a function of n c t, 
where n c and t are the number density of electrons and 
elapsed time, respectively 

The result of the NEIJ model fit is shown in figure 
2-c with the best-fit parameters of kT e \ = 1.07lg Qg kcV, 
fcT e2 = 0.48±0.02keV, and log(n c i [scm" 3 ]) = 11. 83 ±0.05. 
Although the residuals below ~ 2.0 kcV were fairly im- 
proved (x 2 /dof = 928/351), line-like excesses still re- 
mained at ~ 0.8 keV and ~ 1.2 kcV (the arrows in figure 
2-c). 

As pointed out by a number of previous studies on 
SNRs, emissivity prediction of Fe L-shell lines are un- 
certain in currently available plasma code (e.g., Hughes 
et al. 1998; Borkowski et al. 2006; Yamaguchi et al. 2011). 
Gu et al. (2007) claimed that the line intensity ratio 
of FeXVII corresponding to the 3s— >-2p (~ 730 eV) and 
3d— >2p (~ 820 eV) transitions may have a significant er- 
ror. Brickhouse et al. (2000) pointed out that several 
Fe L-shell lines (for example, n > 6^2 for FeXVII and 
n = 6, 7^2 for FeXVIIl) are missing from the SPEX 
code. We thus assumed that the excesses at ~ 0.8 keV 
and ~ 1.2 kcV are caused by such uncertainty, and added 
two Gaussians. Then the fit was improved with a x 2 /dof 
value of 536/349 (figure 2-d). The best-fit parame- 
ters were kT cl = 1.25^;^ keV, kT c2 = 0.48 ± 0.02 kcV, 

log(M [scm~ 3 ]) = 11.85±S;oe- We notc tnat tlie extra 
Gaussians give no significant change of the best-fit param- 
eters of the simple NEIJ model (with no extra Gaussian) . 
Hereafter, we used the NEIJ plus two Gaussians model, 
referring as the NEIJ model for simplicity. 

3.3. Fit with Thermal Plasma plus Power-Law 

We tried the NEIJ model to the other two regions B and 
C. The model gave a nice fit for the region B spectrum, 
but failed to reproduce the region C spectrum with signifi- 
cant excess above -~4.0kcV. This excess is consistent with 

3 More detailed information about this model is found at 
http: / / www.sron.nl/files /HEA/SPEX / manuals / manual.pdf 



the hard band image (figure 1-right; 5.0-7.0 keV), where 
we found a previously unreported feature in the south- 
west part of the FoV. The hard component in the region 
C spectrum was found to be well explained with an ad- 
ditional power-law with a photon index T of 2.12 ±0.11 
(figure 3). 

We finally tried a combined model of the NEIJ and 
power-law components for all the other regions (A and 
B), since we can see a hint of slight hard excess over the 
NEIJ model in the region A spectrum as well (see figure 
2). The power-law index was fixed to the best-fit value 
for the region C (r = 2.12). The best-fit results are shown 
in figure 3 and table 2. We found that the power- low 
component is significant in the region A, while only an 
upper limit was obtained for the region B. 

3.4- Spatial Distribution of the Recombining Plasma 

We further examined more detailed spatial structure of 
the recombining plasma (RP) by dividing the entire FoV 
into 10 small regions and fitting each spectrum with the 
NEIJ plus power-law model. We fixed the photon index of 
the power-law to be 2.12 for all the regions. Statistically 
accepted results were obtained for all the spectra. Figure 
4 shows spatial distributions of the best-fit parameters of 
kT c i (left), kT c2 (middle) and log(n e t) (right). We found 
that the recombining state is achieved in the entire regions 
of the SNR. 

4. Discussion 

4-1. Recombining Plasma 

We found that the spectra from all the regions of W44 
were nicely fitted with the model of the NEIJ and power- 
law components. In order to clarify the property of the 
plasma state, we estimate the average charge per element 
and compared to those expected in a CIE plasma with the 
same electron temperature. The results derived with the 
best-fit parameters for the region A spectrum are shown in 
table 3. We see that the averaged charge states of lighter 
elements, C, N, and O are consistent with those in CIE, 
while those of the heavier elements, such as Si, S and Fe 
are higher than CIE, indicating that these elements are 
in the recombining state. In contrast, some previous X- 
ray studies claimed that the plasma in W44 is ionizing or 
nearly in CIE (Rho et al. 1994; Kawasaki et al. 2005). One 
of the possible origins of this discrepancy is a difference in 
the energy resolution and/or photon statistics in the hard 
X-ray band. 

Our result suggests that either only the initial electron 
temperature kT c \ cooled down rapidly to kT c2 (here, elec- 
tron cooling), or only the initial T z (=kT c i) was raised 
higher than fcT e 2(here, selective ionization). These led the 
ions to be higher charge states (RP) compared to those in 
CIE. In the next sections, we separately discuss these two 
cases. 

4-2. Electron Cooling 

One possible scenario for the electron cooling is ther- 
mal conduction, so we first evaluate this possibility. In 
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Fig. 4. Spatial distributions of the values of kT c \ (left), kT c 2 (middle) and log(n c t) (right). The black contours represent the 
1442.5 MHz radio continuum emission with VLA. The blue contours represent 12 CO (J = 2-1) in the velocity of 40-50 km s _1 
obtained with NANTEN2 4 m submillimetcr telescope of Nagoya University. 



this scenario, the energy transfer from the central region 
(kT e = 0.48 ± 0.01 keV) to the cooler periphery (kT e = 



Table 3. Average Charge Per Element 



0.40+q'q 2 keV) is considered. A scale length of the tem- 
perature gradient Zt = (grad/nT) -1 is calculated to 
be 1.5 x 10 20 cm, at the distance to the SNR of 3kpc. 
Therefore, the thermal conduction timescale i C ond is esti- 
mated to be 



cond 



1 x 10 7 



1.0cm 3 



kT, 



c(in) 



0.48 keV 
h 



1.5 x 10 20 cm 



-5/2 



yr. 



(1) 



(Spitzcr 1962; Kawasaki ct al. 2002), where the electron 
number density n c is calculated from the best-fit EM we 
obtained. The derived timescale is nearly 50 times longer 
than the age of W44 (~ 20,000yr), and hence the thermal 
conduction scenario is difficult to reproduce the observed 
electron temperature gradient. 

The other possibility of the electron cooling is adia- 
batic expansion, which may occur in a rarefaction pro- 
cess; a rapid cooling of electrons occurs when a blast wave 
in a dense circumstcllar medium (CSM) breaks into a 
lower-density ISM (Itoh & Masai 1989; Shimizu et al. 
2012; Moriya 2012). Sawada & Koyama (2012) ar- 
gued the recombination time t Tec is an essential param- 
eter to evaluate the possibility of this scenario. We hence 
estimate the t rcc using the best-fit n c t values: the low- 
est at the east region of 3.8 ±0.1 x 10 11 scm -3 , and the 
average in the center region of ~ 7.1 ±0.1 x 10 1:L scm _3 
(figure 4-c). Since the electron number densities n c at 
the center and east region are, respectively ~ 1.0 cm -3 
and ~ 0.8 cm -3 , the recombining time is derived to be 
~15,000-20,000yr. Wolszczan et al. (1991) estimated the 
age of the pulsar PSRB1853+01, physically associated 
with W44, to be 20,000 yr, which is consistent with our 
estimated recombination time. Therefore, the adiabatic 
cooling scenario is preferable in our case, because the 
shock-heating and the rarefaction should occur at the ini- 
tial phase of the SNR evolution. Broersen et al. (2011) 
derived conditions for general SNRs to keep RP in their 



Element 


Average 
NEIJ (Region A) 


Charge 
CIE at 0.45 keV 


C 


5.99 


5.99 


N 


6.98 


6.98 


O 


7.92 


7.92 


Ne 


9.51 


9.34 


Mg 


10.87 


10.32 


Si 


12.26 


12.01 


S 


13.95 


13.91 


Fe 


16.45 


16.37 



evolutions. Comparing the i roc with an adiabatic cooling 
time t a( j in different n e , they concluded that a low n e sig- 
nificantly decreases the recombination rate. They claimed 
SNR 0506-68 (~ 4,000yr) to be still RP in a dense ISM 
(n c ~ 10cm -3 ). Applying their argument to W44's case, 
we confirmed that the t IOC for W44 is about one order 
of magnitude above that for SNR 0506— 68 (see also Vink 
2012). This estimation also supports the adiabatic cooling 
scenario. 

4-3. Selective Ionization 

Possible mechanisms to rise up kT z (selective ioniza- 
tion) would be either ionization by supra-thermal elec- 
trons (e.g., Tanaka 1986; Masai et al. 2002), or photo- 
ionization (e.g., Piro et al. 2000). In W44, the 2-10 GcV 
and 0.4-3 GeV emissions were discovered with the Fermi 
LAT (Abdo et al. 2010) and AGILE (Giuliani ct al. 2011), 
respectively. The idea of the former scenario is that supra- 
thermal particles are injected into a Fermi acceleration 
process, and are accelerated up to very high energy. The 
supra-thermal particles also ionize atoms in a thermal 
plasma to higher ionized states, and produce the RP. 
In fact, large solar flares often show a signature of RP 
(in Fe atoms) and hard X-ray emissions (Tanaka 1986). 
Thus, the scenario of supra-thermal particle ionization 
may predict positive correlation between the ionization 
temperature (kT z ) and the flux of TeV/GeV 7-rays/hard 
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Fig. 5. 12 CO (J = 2-1) contours of NANTEN2 overlaid with 
the hard X-ray image in the 5.0-7.0 keV band. The blue and 
yellow contours are the 12 CO fluxes in the velocity ranges 
40-50 km s — 1 and 50-60 km s -1 , respectively. The white 
square is the XIS FoV. 

X-rays. However, our observation shows no systematic 
spatial variation of kT z , and hence possibility of selective 
ionization by supra-thermal particles is unlikely. 

Another possibility, photo-ionization, requires a strong 
X-ray source in or near W44, but none of such sources is 
found at present. Therefore, such ionizing source should, 
if any, be bright in past. A 7-ray burst is one of the 
candidates. If W44 was a remnant of hyper- nova explosion 
associated with a 7-ray burst, the X-ray afterglow would 
be strong enough to ionize a substantial amount of atoms 
(high kT z — kT e i), whereas the electron temperature was 
relatively low (low kT C 2). A weak point of this scenario 
is, however, that no particular signature for a hyper-novae 
has so far been observed in W44. 

4-. 4- Origin of Hard X-ray Emission 

We detected strong hard X-rays from the west re- 
gion of W44 (region C), and possibly from the entire 
SNR regions. A conceivable origin is a pulsar wind neb- 
ula (PWN) associated with PSRB1853+01 (Wolszczan 
et al. 1991; Petre et al. 2002) located just outside the 
FoV (the white cross in figure 1). The Chandra ACIS 
found that the X-ray emitting region of the PWN is 
within - 2' from PSRB1853+01, while we found that 
the peak of the hard X-rays is far from the PWN (~ 
12'). The total flux of the PWN plus PSRB1853+01 is 
~ 1.7 x 10~ 13 ergscm -2 s _1 in the 2-8 keV band (Petre 
et al. 2002), far less than the hard X-ray flux in the region 
C (~ 1.37 x 10~ 12 ergscm~ 2 s _1 ). We also simulated pos- 
sible contamination of the strayed flux from the PWN and 
PSRB1853+01 to the region A by using xissim (Ishisaki 
et al. 2007), and confirmed it is negligible. We thus con- 
clude that the hard X-ray emission we discovered is inde- 
pendent from the PWN or PSRB1853+01. 



The hard X-rays was found to have an arc-like structure 
(figure 1-right) and are correlated with a filament of the 
radio continuum (1442.5MHz with VLA). For more de- 
tailed study, we compared the spatial distributions of the 
hard X-rays with ambient molecular clouds, 12 CO (J = 2- 
1) emission obtained by NANTEN2. The result is shown 
in figure 5. We see a clear anti-correlation between the 
hard X-rays and the CO cloud in the velocity range 50- 
60kms -1 (the yellow contours). This cloud is identical 
to what was observed by Nobeyama 45-m radio telescope, 
Cll named by Seta et al. (2004), and was suggested to 
be interacting with the SNRs blast wave (see figure 11 in 
their paper). 

The arc-like structure of the hard X-rays suggests that 
this emission comes from the shock. A typical shock veloc- 
ity of W44 is - 650 (fcT G /0.5keV) - 5 kms" 1 . As far as we 
assume a simple direct diffusive shock acceleration, such 
slow shock is insufficient to explain the flux of the observed 
hard X-rays. We presume that the flux enhancement pos- 
sibly arises in association with the shock-cloud interaction. 
We focus on pre-existing cosmic ray (CR) electrons com- 
pressed and reaccelerated between the shock and dense 
cloud. In the radio band, the resultant synchrotron en- 
hancement is known as the van der Laan mechanism (van 
de Laan 1962a; 1962b) which accounts for the radiative 
radio shell (Blandford & Cowie 1982). Based on a multi- 
wavelength observation of W44, Cox et al. (1999) pointed 
out that the radio synchrotron emission is possibly en- 
hanced by the pre-existing CR electrons which are reac- 
celerated by the van der Laan mechanism. Uchiyama et al. 
(2010) applied a similar mechanism to some middle-aged 
SNRs including W44 to interpret their GeV 7-ray emis- 
sion. They concluded that it is sufficient for the reacceler- 
ating of the pre-existing CR through the hadronic process 
to explain the observed GeV 7-ray intensity. We simply 
expanded their argument to explain the observed hard X- 
rays. By our rough simulation, if the hard X-ray is derived 
from such mechanism, the electrons shall be accelerated 
to at least ~100TeV. In this case, the acceleration time 
ta.ee is calculated to be 

VIOOTcVAloAlOOkms- 1 / 

where cp m ax is the maximum attainable energy, 77 is the 
gyrofactor (Uchiyama et al. 2007), v s and Bq are the shock 
velocity decelerated by cloud and the pre-shock magnetic 
field (Reach et al. 2005), respectively. These parameters 
were set as Uchiyama et al. (2010) had assumed for W44 
in their calculation. As a result, the t acc becomes signifi- 
cantly longer than the age. 

A local anti-correlation between CO clouds and non- 
thermal X-ray emissions is also found in another SNR, 
RXJ1713. 7-3946 (Sano et al. 2010). Performing magne- 
tohydrodynamic (MHD) simulations of a SNR shock pass- 
ing through clumpy interstellar clouds, Inoue et al. (2012) 
successfully interpreted this result; the synchrotron X-ray 
enhancement occurs in spatially small regions at the vicin- 
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ity of the cloud due to an amplified magnetic field caused 
by turbulent dynamo action (see also, Giacalone & Jokipii 
2007). One major difference between RX J1713. 7-3946 
and W44 is a shock velocity. Inoue et al. (2010) showed 
that a slow shock wave (~ 500 km s _1 ) attains a maximum 
magnetic field strength of ~ 400 /J.G. Given a typical shock 
velocity of W44, ~ 650 km s , however, the resultant i a cc 
becomes substantially longer than the synchrotron cooling 
time t syn ; t syn /t acc < 0.07. We concluded that this scenario 
is also difficult to explain the observed X-ray flux. 

In either case, the acceleration time t acc is too long to 
obtain TeV electrons. Alternatively, if TeV protons are 
still sufficiently remaining in W44, the hadronic process 
provides secondary electrons which may peak at TeV en- 
ergies. As already claimed by Abdo et al. (2010) and 
Giuliani et al. (2011), the GeV 7-ray in W44 is possibly 
due to the hadronic process. We found no spatial cor- 
relation between the hard X-ray emitting region and the 
global GeV 7-ray distribution (Abdo et al. 2010; Giuliani 
et al. 2011). The reason would be that the synchrotron 
X-ray is enhanced only around a shocked clumpy cloud 
where the magnetic field is amplified by the mechanism of 
Inoue et al. (2010). The remaining problem is a CR escape 
which determines a peak energy of the existing protons, 
hence the secondary electrons in W44. While the max- 
imum possible energy differs by more than one order of 
magnitude among theories (e.g., Ptsukin & Zirakashvili 
2003), it is generally difficult for a middle-aged SNR to 
have 100- TeV order of protons. The origin of the hard 
X-rays remains an open question. 

5. Summary 

1. The thermal X-ray emission in W44 is centrally- 
filled, similar to other MM SNRs, while the arc- 
like hard X-ray emission has been found for the first 
time. 

2. The thermal X-ray spectra are well represented with 
a recombining plasma (NEIJ model). 

3. The electron temperature of the thermal plasma 
would be cooled down rapidly from ~ 1 keV to 0.4- 
0.5 kcV at -20,000 years ago. 

4. The most plausible origin of the recombining plasma 
is rarefaction in the early phase of the SNR evolu- 
tion. 

5. The surface brightness of the hard X-ray and 12 CO 
(J = 2-1) emissions are spatially anti-corrclated. 

6. Some particle acceleration mechanisms which we 
discussed give no plausible explanation of the anti- 
correlation. The origin of the hard X-ray emission 
is an open question. 
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